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Light-like galileon solutions have been used to investigate the chronology problem in galileon-like 
theories, and in some cases may also be considered as solitons, evading a non-existence constraint 
from a zero-mode argument. Their stabilities have been analyzed via "local" approximation, which 
appears to suggest that all these light-like solutions are stable. We re-analyze the stability problem 
by solving the linear perturbation equation exactly, and point out that the finite energy condition 
is essential for the light-like solitons to be stable. We also clarify potential ghost instabilities and 
why the zero-mode argument can not be naively generalized to include the light-like solitons. 



I. INTRODUCTION 



The galileon model is a scalar field theory that pos- 
sesses a novel field symmetry tt — > tt + 6^2;'^ + c, fe^, c = 
const [ij . It is originally constructed to describe large dis- 
tance modifications of General Relativity 1], which can 
capture some local 4D gravitational physics of certain 
braneworld scenarios 043] as well as the non-linear mas- 
sive gravity model @ (see 0, Q and references therein 
for details and related cosmology). Albeit simple, this 
novel field theory had been rarely explored before, but 
has been since found to have interesting field-theoretic 
properties. The general galileon Lagrangian in (n + 1)-D 
is given by 



n+l 



_ (1) 

where [ ] means anti-symmetrization of the indices and 
Qfi's are the couplings of the theory. Despite the ap- 
parent high order derivatives, its equation of motion re- 
mains second order, avoiding the Ostrogradski ghosts. 
Also, perhaps more interestingly, all the couplings in this 
model are not renormalized by loop corrections [1, Q. 
So it is interesting to investigate its classical non-linear 
dynamics, particularly to search for stable solitonic so- 
lutions. As the non-linear galileon terms contain higher 
order gradients, it may be possible that Derrick's the- 
orem [lOj can be evaded. However, using a zero-mode 
argument, Endlich et al fllj have shown that there is no 
stable static soliton in the single galileon model. This 
may be evaded by a multi-field generalization coupled 
with a non- linear sigma model constraint [T2| . along the 
lines of the skyrmions [l^ [13] , which have been used to 
model baryons and resonance states. 

Another way to evade the zero-mode constraint, as 
recently pointed out by Masoumi and Xiao is to 

construct solutions that can not be boosted to a rest 
frame, i.e., a solution that is light-like or wave-like (see 



also |16l4l8l |). In 1+lD ^, the only non-linear galileon 
term is the a2 coupling, in which case it has been found 
that any left or right moving mode tt = ttq (a; ± t) is a 
solution to the equation of motion (TB - fTsl . The reason 
for this can be understood via the "target space diffeo- 
morphism" symmetry fl9, 20]. Or, one may simply see it 
from the equation of motion in the light cone coordinates 

{u — X + t,V = X — t), 



dudvTT + 2a [{dud^n)'^ - S^ttS^tt] = 0, 



(2) 



where a = 3a2/4. The null modes n — t:q{x ± t) vanish 
the linear and non-linear pieces simultaneously. If the 
solution (the field derivative, not necessary the field it- 
self, as will be clear later) is a concentrated wavepacket, 
we may consider it as a soliton. This is in contrast to a 
typical soliton where the non-linear terms cancel the lin- 
ear terms in the equation of motion. Therefore, galileon 
solitons can form even when the non-linear interactions 
are small. 

These light-like solutions have been used to investi- 
gate closed time-like curves and their implications for the 
chronology protection problem in various theories [1^ 
[T8| . However, so far, the stabilities of these light-like so- 
lutions have only been tackled in the "local" approxima- 
tion, i.e., by calculating the eigenvalues of the kinetic co- 
efficient matrix of the perturbative Lagrangian, based on 
which the light-like solutions appear to be stable for any 
light-like configurations n = no{x ± t) [l^ [l^, putting 
them on the same footing as the linear wave equation. 
As one may expect, this "local" stability analysis is by 
no means rigorous (see e.g. [1, (U), and this peculiar 
conclusion may be merely an artifact of this approxi- 
mation. This, however, can only be determined if the 
linear perturbation equation in question can be solved 
exactly, which typically is difficult. In this note, we will 
re-examine the stability problem of the light-like galileon 
solutions, and show that the relevant linear perturbation 
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^ In this note, we will mainly focus on 1+lD for simplicity. In 
3-l-lD one may equivalently consider homogeneous 2-brane or 
domain wall configurations. 
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equation is fortunately solvable. The finite energy con- 
dition of tlie given solution, which has been largely ne- 
glected so far, is important for the solution to be stable. 
There may be tachyonic instabilities (or growing modes) 
if the light-like configurations are not properly chosen to 
have finite energies. Also, we will point out that ghost 
instabilities will develop for certain configurations, which 
has been neglected in [15| , and clarify why the zero- mode 
argument of [11] can not be naively generalized to include 
the light-like soliton case. 



II. LINEAR STABILITY OF THE LIGHT-LIKE 
SOLUTIONS 

As discussed in the Introduction, any light-like config- 
uration Tro{x zt t) is a solution to the galileon equation 
of motion Since we are only interested in the linear 
stability, we perturb a given solution t: — ttq^x ± t) + (jj 
and obtain the pcrturbative action 



(3) 



The kinetic coefficients are 



1 + a{x ± t) ^a{x ± t) 
^a{x±t) -l + a{x±t) 



(4) 



where we have defined the function a{x ± t) = 2aTTo{x ± 
t). From the pcrturbative action ([3]), we get the linear 
pcrturbative equation of motion: 



(5) 



where we have used the relation d^Z^^'^ — 0, thanks to the 
galileon construction. So the problem is to solve Eq. ([5]) 
exactly. Note that Eq. ^ is of the hyperbolic type, since 
the discriminant A = - {Z^'^f = -1 < 0, so it 

has a well-defined Cauchy initial problem. 



A. Tachyonic instabilities 

First, we will be looking at the tachyonic instabili- 
ties, that is, whether there are modes that grows (quasi- 
) exponentially with time. Unable to solve the second 
order partial differential equation exactly, this is often 
tackled in the "local" approximation, i.e., by assuming 
the background varies slowly both spatially and tempo- 
rally so that the kinetic coefficients Z^^'^ can be treated as 
constants. This is justified if we are only concerned about 
the ultraviolet modes: When we "zoom in" at a given 
spacetime point, the smooth functions Z'"^(a;'^) certainly 
vary slowly. Then the linear stability problem is reduced 
to determining whether Z^'^ has the right Lorentz signa- 
ture at every spacetime point, or, equivalently whether 
^•oj^-oi _ z^^Z''^ is positive definite For the problem 



at hand, Z^^'^ indeed has the right signature, independent 
of the form oi a{x ±t) = 2a7fb(a; ± <), as observed by 
[l^[2l|. However, this "local" approximation may miss 
some important instabilities such as resonances with the 
background. For the problem at hand, as one may ex- 
pect, the stability problem being independent of a{x ± t) 
could well only be an artifact of this approximation. So 
solving the linear perturbation problem exactly^ if possi- 
ble, is certainly desirable, and fortunately, is achievable 
for the problem in question. 

We will solve it in the light cone coordinates (u = x + t, 
V = x — t), where the pcrturbative equation of motion can 
be written as. 



a{u)d^(f> - dudy(f) = 
a{v)dl(j) - dudy<j> = 



for ttq = nQ{x + t), (6) 
for TTo = 7ro(a; - i), (7) 



respectively for the left and right moving modes. As 
expected, these equations are symmetric in terms of ex- 
changing u and V. 

Without loss of generality, we can focus on the right 
moving mode ttq — '7To{x — t), that is, Eq. ([7]). To solve 
Eq. ([7]), it is convenient to extend (j) to be complex valued 
and separate the variables (l){u,v) = U{u)V{v). This 
leads to 



dlU{u) d,V{v) 



duUiu) a{v)V{v) 



A, 



(8) 



where A is a constant, or an eigenvalue for the relevant 
eigenmode. Note that here we should allow A to sample 
the whole complex plane. Now, we can solve U{u) and 
V{v) separately, whose general solutions are respectively 
given by 



U{u) = Ci + 026^ 



(9) 
(10) 



where ci, C2, C3 are integration constants and the function 
b{v) = / a{v')dv' = 2aa„7ro(w). (11) 



Thus, in the (i, x) coordinates, the eigenmodes of (j) are 

(12) 

where 

B{x -t) = b{x -t) + x - t. (13) 

At this stage, Ai and A2 are arbitrary complex numbers 
and the most general solution to Eq. ([7]) is the linear 
combination of all the these eigenmodes with different 
Ai and A2, which can be written as 



(l>it,x) 



dX 



Pi(A)0i')(i,x)+p2(A)0f (i,a;)l , (14) 



where pi and p2 are distribution functions of the eigen- 
modes, yet to be constrained by the initial and boundary 
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conditions of tfie given problem. Only eigenmodes that 
can pass the "filtering" of the initial and boundary con- 
ditions can be used to construct the solution space of the 
given problem. If there are (quasi-)exponentially growing 
modes that survive the "filtering" , then the background 
solution in question is unstable. Otherwise, it is free of 
tachyonic instabilities. 

The initial and boundary conditions in question are 
that the perturbation field (f){t, x) and its "velocity" 
(j){t, x) should be small, or bounded, initially. Since nei- 
ther the eigenmodes of (j){t, x) nor those of (f>{t, x) de- 
generate at t = 0, this requires that every eigenmode 
of ^f^{t,x) and its "velocity" ^f^it, x) are bounded at 
t = 0. (Note that degeneration of eigenmodes at i = 
happens, for example, for the case of the 1-l-lD wave 
equation, in which case the reduction to the eigenmodes 
is not quite appropriate. We will discuss this in the Ap- 
pendix.). Further, since 

0^ ^ 0, x) and {t = 0, x) 
are smooth functions, it is enough to require (j)"^, {t = 

0, a; — > ±cxd) and <t>{^.{t = 0,a; — >■ ±oo) to be bounded, 
and these are our "filters" . Any eigenmode, either itself 
or its "velocity" , that is not bounded initially should be 
thrown away. The eigenmodes and their "velocities" at 
time t — are 

0W(O,x)=e^^H.), 0f^)(O,.) = e^^[^(-)l, (15) 

0^')(O,x) oc a(a;)e^i*("\ 0i'^(O,x) oc (a(x) - l)e^^[^(")l, 

(16) 

where, as before, a(x) = h'(x) and B(x) — b{x) + x. So 
whether these eigenmodes can pass through the "filters" , 
consequently the stability of the problem, depends on the 
form of the background configuration 6(a;), as well as the 
values of Ai and A2. 

Note that if there are eigenmodes with pure imaginary 
Ai and A2 passing though the "filters", then the per- 
turbations will oscillate around the background solution, 

1. e., the given solution will slightly "breath in and out" 
once perturbed. This is like a "limit cycle" in a standard 
dynamical system, which is regarded as stable. Potential 
instabilities arise only when Ai and A2 with non-zero real 
parts can pass though the "filters" . 

Let us first focus on the first "filter": the eigenmodes 
4>^^.{t,x) should be bounded at (i = 0,a; — J- ±00). 

First, we note that whenever 0^^^^(O,a;) — e^^''^^^ is 
bounded, (j)"^'^{t^x) — e'^i''^^^*) is also bounded, so there 

are no growing modes from the 0^^^ {t, x) sector. 

(2) 

The {t, x) sector is more complicated. Since the 
first "filter" requires e^''^^''\ hence also e^2-B(x-t)^ 
bounded, we see from Eq. that modes with Re{X2) > 
are (quasi-) exponentially growing. (The potential loop- 
hole for this argument is that, even for i?e(A2) > 0, when 
t — ?• -l-oo, it is possible that e^2B(x-t) _i. q at a rate such 
that ^'^^^'tQ^^B(x-t) bounded. However, this is only 
true for certain regions of x. Indeed, as B{x — <) is a 



regular function oi x — t, there are always regions of x 
such that B{x — t), hence also ^^-^Bix-t) ^ remains finite.) 
If only modes with Re{\2) < can pass through the first 
"filter" , without the need of examining the second "fil- 
ter" , we can say that there are no tachyonic modes for 
the given configuration b{x — t). To proceed, let us clas- 
sify the initial configuration B{x) according to different 
asymptotes: 

1. B{x — s> ±cxd) — ?> —00, or one end is bounded and 
the other goes to ~oo: modes with Re{\2) > can 
pass though, thus this case is potentially unstable. 

2. B{x — > ±cxd) — > +00, or one end is bounded and the 
other goes to +00: only modes with i?e(A2) < 
can pass though, thus this case is stable, without 
the need of examining the second "filter" . 

3. B{x — J- ±00) — > =Fc»- only modes with i?e(A2) = 
can pass though, thus this case is stable, without 
the need of examining the second "filter" . 

4. B{x — >■ ±00) are both bounded: modes with any 
A2 can pass though, thus this case is potentially 
unstable. 



Now, let us consider the second 



"filter": ^f^{t,x) 
should be bounded at (t = 0,a: — > ±00). From Eq. (|16l) . 
we can see that this "filter" is equivalent to requiring that 



dx 



±00] 



and dx 



0i'j(O,a; = ±oo) (17) 



are bounded. Note that this is not implied by the re- 
quirement that (0, X = ±00) and (j)"^^ (0, x = ±00) 
are bounded, as one would have naively thought. For 
example, while sin(x^) is bounded within {—00, -l-c»), its 
derivative is not. Further classification of B{x) is rather 
complicated, whose details we shall not go into. We just 
point out that, depending on the background solution 
b{x), the second "filter" may kill off more unstable modes, 
as well as stable modes, which have passed through the 
first "filter". Once the form of the background solution 
is given, the condition (jl7p can be readily checked and 
the full stability problem is therefore easily determined. 

So we see that not all light-like (or wave-like) solutions 
■n ~ ■n(){x — t) are stable, in contrast to what one would 
infer from the analysis of the "local" approximation. 

So far, we have been considering general solutions, so- 
lutions that may or may not have non-trivial global struc- 
tures. This is important if one is interested in, for exam- 
ple, globally oscillatory solutions or processes involving 
extreme long wavelengths. However, solitons are nor- 
mally defined as solutions that are localized and have 
finite energies. If we are interested in localized solutions 
with finite energies, then we will see that galileon solitons 
in 1-l-lD are always stable. Note that it is not necessarily 
for the galileon field itself to be localized, but its deriva- 
tive has to fall off quickly enough so that the energy 



E[M = 2 j Av{d^TTo{v)f ^^Jdv b{v 



(18) 
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is finite. This requires b{v) to fall off quickly at infinities, 
so we have B(x — >■ ±00) — x. Therefore, the finite energy 
soliton case corresponds to the case 3 above, which is 
fortunately stable. From the first "filter', we know that 
A2 has to be pure imaginary in this case. Also, since b{v) 
has to fall off quickly at infinities, the second "filter" is 
automatically satisfied. So, if perturbed, galileon solitons 
will "breath in and out" during their evolutions. 

B. Ghost instabilities 

The kinetic term of the action (O is = (1 + 

a{v))(jP (Again, without lost of generality, we have fo- 
cused on the right moving mode.). Clearly, when a{v) < 
— 1, there is a ghost instability, which was neglected in 
[l5| . To avoid the ghost instability, we require 

a{v) = 2a7fo(i') > -1. (19) 

This does not contradict the fact that Z'^'^ is negative def- 
inite, because the signatures of the temporal and spatial 
components can swap. If the scalar tt is a standing only 
field, a ghost instability is not really a problem per se. 
After all, the action can always be re-defined by multiply- 
ing an overall sign, which does not change the equation of 
motion and leads to a well-defined path integral formula- 
tion. But when tt couples to other fields such as gravity, 
a ghost signals very serious instabilities, particularly in 
the quantum theory (see e.g. (2^). 

III. EVASION OF THE ZERO-MODE 
ARGUMENT 

In [15| , it was pointed out that the zero-mode argument 
of [ll| does not include the case where the soliton can 
not be brought to a rest frame. But one might wonder 
whether we can apply a similar instability argument to 
the light-like solitons. After all, a soliton moving at the 
speed of light can still be invariantly translated spatially 
or temporally. Here we explicitly show why similar zero 
mode arguments can not be naively applied to the light- 
like solitons. 

First, we review Endlich et aFs zero-mode argument for 
static solitons in the single galileon theory [11]. Suppose 
there exists a solitonic solution which can be boosted to 
be static tt = ttq{x^) (In this section we work on general 
n 4- ID dimensions.). To decide whether this soliton is 
stable or not, one perturbs the soliton tt = 7ro(a;*) -I- 
and see whether the energy of the soliton can be lowered 
by some perturbation (p. For the galileon theory, the 
perturbative energy of the static soliton is given by 

52H = J d''x{(})^ + Z'^d^(j)dj(b), (20) 

where 

Z^3 = _ 2a{dkd''TToS'^ - d'd^TTo) + ■■■ . (21) 



Since the energy of the soliton is invariant under spatial 
translations, the perturbation 

= TToix' + e') - TToixl = e'9,7ro + 0{e^) (22) 

should give rise to = Q + 0{e'^), i.e., this is a zero en- 
ergy mode. Now, we can use this zero mode to probe the 
positivity of the kinetic coefficients Z^^ . Far away from 
the center of the soliton, the definition of solitons requires 
Ottq, ddiTo to be very small, so we have Z^^ ~ 6^^ . This 
implies that S2H gets positive O(e^) contribution from 
the integration far away from the center of the soliton. 
For the zero energy mode, there must be some negative 
contribution near the center of the soliton, which im- 
plies that Z'-' has negative eigenvalues near the center. 
The perturbation 0, hence di4>, is essentially initial and 
boundary conditions of the system that can be arbitrary 
chosen. If di4> are chosen to be deformations of the zero 
mode that "enhance" the negative eigenvalue directions 
of Z'^^ , we can get negative S2H. Therefore, the assumed 
soliton is not stable ^ . 

Now, why can we not develop a similar argument for 
the null solitons, using the spacetime translation sym- 
metry the system still possesses? For the static soliton 
case, we were actually a bit sloppy in differentiating the 
perturbative Hamiltonian (Hamiltonian derived from the 
quadratic perturbative Lagrangian) and the full Hamil- 
tonian (Hamiltonian derived from the full Lagrangian 
and then perturbed to quadratic order). The space- 
time translation invariance is a symmetry of Minkowski 
space, i.e., a symmetry of the full Lagrangian or Hamil- 
tonian, but not necessarily a symmetry for the pertur- 
bative Lagrangian or Hamiltonian. So the full Hamil- 
tonian is the one we should really be using, if we want 
to make use of the spacetime translation symmetry. On 
the other hand, the zero-mode argument above relics on 
the use of the quadratic Hamiltonian. This is well jus- 
tified as long as the background solution, which mini- 
mizes the Lagrangian, also minimizes the Hamiltonian, in 
which case the quadratic perturbative Hamiltonian and 
the full Hamiltonian only differ by a constant, the en- 
ergy of the background solution. Whilst this is the case 
for a static background, it is not necessarily true for a 
time-dependent one such as a light-like soliton. There- 
fore, for a time-dependent background, apart from the 
zero-order piece and the quadratic piece, the full Hamil- 
tonian may include a linear piece, and the zero-mode as- 
sociated with the space translation does not necessarily 
vanish the quadratic Hamiltonian (up to second order in 
perturbation). Instead, it should vanish the linear and 
quadratic pieces combined. Therefore, the zero-mode ar- 
gument generally does not apply to a time-dependent 
background. 



^ In 1 + ID, there is actually no non-trivial static solutions, as the 
equation of motion reduces to a trivial linear equation. 
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We can check that this is indeed the case for our situa- 
tion at hand. The full Hamiltonian of the I+ID galileon 
theory is given by 



H = J dx[{l + 2a7r")7r^ + vr'' 



(23) 



Perturbing it (tt = 7ro(a; — t) + 4>(t,x)) up to quadratic 
order and integrating by parts, we have 



H — Hq + Hi + H2 



(24) 



with 



Ho = J da;27fo^ 

Hi= dx[2{l + 2a7fo)7fo0 - 2(1 - 2a7fo)7fo^'], 



H2 = I dx[{l + 2ano)(j)^ + (j>'^ + Aano^cj)"] 



Substituting in the zero-mode (j) — T^oix — t + e) — 
■7Ta{x — t) = — eTTo -|- e^7fb/2 -|- 0{e^), integrating by parts 
and collecting the 0{e^) terms, we have Hi + 0{e^) = 
-2e2/dx7fo' = -i?2 + 0(e')- 



IV. DISCUSSION 

In this note we have re-examined the (linear) stability 
problem of the light-like galileon solutions. These solu- 
tions may be considered as solitons in some cases and are 
interesting in their own right. They have also been used 
to investigate closed time-like curves and their implica- 
tions for the chronology protection problem in various 
galileon- like modified gravity theories [Tsl - fisl . By solving 
the second order linear perturbation equation explicitly, 
we have been able to determine whether there are (quasi- 
) exponentially growing modes for given background con- 
figurations and classify general configurations according 
to their behaviors at the initial infinite boundaries. For 
the case of a finite energy soliton (or wavepacket), it is 
the finite energy condition of the soliton that guaran- 
tees the stability of the soliton, a point that has been 
overlooked so far. For a general solution with non-trivial 
global structures, the stability of the solution is more 
subtle. But if the configuration is given, its stability can 
be readily checked using the conditions established. Ex- 
amining the positivity of the kinetic coefficient matrix 
is frequently used to determine the absence of tachyonic 
instabilities in linear perturbation theory. This mathe- 
matically solvable case exemplifies the expected limita- 
tion of this approach. We have also pointed out that 
for certain solitonic configurations ghost instabilities can 
also develop, which will become a serious problem when 
the galileon is coupled to other fields, and shown why 
the zero-mode argument of Endlich et al |Tl| can not be 
naively generalized to include the light-like solitons. 



Finally, we pass by a couple of curious observations 
about the galileon solitons. In this simple non-linear the- 
ory, we have a situation that localized, or wavepacket, 
solutions are stable while very long wavelength solutions 
are often susceptive to disintegration. If a system starts 
with a long wavelength dominated configuration, small 
wavepackets, or "massless particles", can emerge from 
it during its evolution, which may have interesting late 
time chaotic phenomena. This is in contrast to the lin- 
ear wave equation, in which case even globally non-trivial 
solutions are stable. Also, the theory with only the a2 
coupling is exactly the decoupling limit of the braneworld 
DGP model 0, Q, where the extrinsic curvature of the 
DGP brane is roughly K^j_^ ~ d^d^ir. For the light-like 
solution, we have K = ry'^^if^^ ~ □7ro(a; ±t) = 0, whose 
geometric origin may be interesting to investigate. 



Appendix A: Stabilities of the 1+lD Wave Equation 

Any function "00(2^ ± t) is a solution to the 1-l-lD wave 
equation tp — ip" = 0. The perturbation on this solution 
(p — -ip — is also governed by the wave equation 

'f-ip" = 0. (Al) 

Its eigenmodes are ip^^'^{t,x) — e'*'i(^+*) and ip^^^{t,x) — 
gA2(a:-t)^ so the general solution is given by 



ip{t,x) 



dA 



/i(A)(pi'^-f/2(A)(pf 



(A2) 



To decide its stability, we can let this general solution to 
pass through the "filtering" of the initial and boundary 
conditions. 

We first focus on the condition that </?(* = 0, x = ±cx)) 

is bounded. Now, the eigenmodes ip^^'^ {t, x) and ip'f'' (t, x) 
degenerate at t = 0. In this case, we may not simply re- 
duce the initial and boundary conditions for Eq. (IA2p to 
a bunch of initial and boundary conditions for the eigen- 
modes individually. To see this, we do the "filtering" 
in the two ways and check whether they are consistent. 
First, impose the initial and boundary conditions that 
the eigenmodes (p^^^ (0, x = ±00) = ip^^^ (0, x = ±00) are 
bounded. This essentially restricts A to be pure imag- 
inary, and therefore it appears that all unstable eigen- 
modes have been "filtered" away. On the other hand, 
if we only require that (p{t = 0,x = ±00) are bounded, 
there are unstable eigenmodes that can survive this "fil- 
tering" . This can be achieved by designing special mode 
distributions. For example, if /i(A) = —/2(A) is satis- 
fied within (— c, +c), c being a positive real number, then 
modes with A € (— c, +c) can also survive the "filtering" 
ofip{t — 0,x — ±00) being bounded. But, from Eq. (jA2[) . 
these modes are clearly unstable in terms of time evolu- 
tion. Therefore, we see the two ways of "filtering" are 
inequivalent. 

Nonetheless, any solution ipo{x ± t) to the wave equa- 
tion is (by accident) actually stable, which can be seen 
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by taking into account the second "filtering" condition 
that the field's initial "velocity" is also bounded. 
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